Error bounds in MEG multipole localization 
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1 Introduction 

Magnetoencephalography (MEG) is a non-inva- 
sive method that enables the measurement of the 
magnetic field produced by neural current sources 
within the human brain. Unfortunately, MEG source 
estimation is a severely ill-posed inverse problem. 
The two major approaches used to tackle this prob¬ 
lem are “imaging” and “model-based” methods. The 
first class of methods relies on a tessellation of the 
cortex, assigning an elemental current source to each 
area element and solving the linear inverse problem. 
Accurate tessellations lead to a highly underdeter¬ 
mined problem, and regularized linear methods lead 
to very smooth current distributions. An alternative 
approach widely used is a parametric representation 
of the neural source. Such model-based methods 
include the classic equivalent current dipole (ECD) 
and its multiple current dipole extension [1]. The def¬ 
inition of such models has been based on the assump¬ 
tion that the underlying sources are focal and small in 
number. 

An alternative approach reviewed in [4], [5] is to 
extend the parametric source representations within 
the model-based framework to allow for distributed 
sources. The multipolar expansion of the magnetic 
field about the centroid of a distributed source readily 
offers an elegant parametric model, which collapses 
to a dipole model in the limiting case and includes 
higher order terms in the case of a spatially extended 
source. While multipolar expansions have been 
applied to magnetocardiography (MCG) source 
modeling [2], their use in MEG has been restricted to 
simplified models [7]. The physiological interpreta¬ 
tion of these higher-order components in non-intui- 
tive, therefore limiting their application in this 
community (cf. [8]). 

In this study we investigate both the applicability 
of dipolar and multipolar models to cortical patches, 
and the accuracy with which we can locate these 
sources. We use a combination of Monte Carlo anal¬ 
yses and Cramer-Rao lower bounds (CRLBs), paral¬ 
leling the work in [3] for the ECD. Results are 
presented for both point sources and cortical patches. 
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2 Methods 

2.1 MEG multipolar expansions 

The external magnetic field is generated by the 
sum of the primary neural activity , designated by the 
current density vector f (#•'), and the volume or return 
currents resulting from the electric field produced by 
the current source. It is the primary currents that are 
the sources of interest in MEG inverse problems [1]. 
In the special case treated here of radial measure¬ 
ments for sources confined to a spherical volume, the 
volume currents do not contribute to the measured 
field, and the radial component bfr) of the magnetic 
field b(r) at location r is given by the well known 
equation: 

b(r) = = r . . ^ f M(r')/cf(r, r')dr' (1) 

r r r 4n*v 

where d(r , r') = r- r' is the distance vector between 
the sensor and source locations, d(r, r') = ||r-r'|| the 
corresponding scalar distance, r/r is a unit vector 
pointing in the radial direction, and V is any volume 
containing the primary source activity. We define the 
magnetic moment density or magnetization as 
M(r') = r' x/(r'). 

The first order multipolar representation is 
derived using a truncated Taylor series expansion of 
the distance d(r,r ') about r /9 the centroid of the 
region to which the primary source is confined: 

d(r , r l + x) 3 = d(r, rj) 3 + 3 d(r, rj) 5 (x • d(r , r ()). (2) 
We truncate the series after the lst-order term, with 
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the assumption that x « d(r, r() , which is to say that 
the source does not have a large radius relative to the 
distance to the observation point. Inserting (2) into 
(1) yields the magnetic multipolar expansion 
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where the integration is carried out over the volume 
of the primary source activity, centered on r { . 
Equivalent Magnetic Moment (EMM): 

Let the extent of the primary source activity V be 
sufficiently small that the second term in the multipo- 



lar expansion is negligible. As shown in [4], we can 
rewrite (3) as 
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where we define m to be the magnetic moment 


m{rj) = j(r l + x)xf(r, + x)dx. (5) 

Furthermore, if we define q{r t ) = jf(r/ + x)dx to be 
the equivalent current dipole , we can express (5) as 
m (r t ) = r t x #(r^) + ^(r,) (6) 


where ih(rj) = Jjc x/*(r z + jt)<A is the magnetic dipole 
moment centered at r, due to the primary current 
density. Thus this model includes the equivalent cur¬ 
rent dipole as a limiting case where the source either 
has virtually no spatial extent or no net magnetic 
dipole moment. 

First-Order Multipole: 

If the spatial extent of the source is sufficiently 
large, then we retain the first two terms in the Taylor 
series and rewrite (3) as 

, KQini-d^n))) _ 
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where Q(rj) is the magnetic quadrupolar term defined 
as the 3 x 3 tensor product Qirj) = [M(r / + x)xdx . We 
can rewrite this tensor using the Kronecker product 
a® b , defined as the concatenation of the product of 
each element of a with the vector b , and the operator 
vec(^), defined as the concatenation of the columns 
of a matrix A into a vector: 

< r 3 (d(r,r,)®r) ) 

+ -- • vec(0 (8) 
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where for notational convenience we drop the depen¬ 
dency of m and Q on r r We therefore characterize 
the first-order multipole using the combination of the 
three magnetic moment vector m , the nine magnetic 
quadrupolar terms in Q , and the location r t . 


2.2 Cramer-Rao lower error bounds 


The Cramer-Rao lower bound (CRLB) is an 
important result in estimation theory that establishes 
a lower bound on the variance of any unbiased esti¬ 
mator of a set of unknown parameters. We model the 
MEG data as 

F=G(l)u + N (9) 

where G(l) is the gain matrix at location /, and u is 
the magnitude and orientation of the source moments, 
i.e. u = q or u = m for the ECD and EMM models 
respectively. We restricted our theoretical error anal¬ 
ysis in this paper to the dipolar models (see comments 
in Section 4. iV is the additive noise assumed to be 
white and Gaussian. We denote the parameters to be 
estimated by 6 = {lu} . The CRLB is given by the 
inverse of the Fisher Information Matrix (FIM): 


J,, = ICRLBmy' 


d(Gu) 

ae. 


( 10 ) 


for i,j= 1, ..., P ; where P is the number of parameters 
and v is the noise variance. 


In this study, we are specifically interested in the 
additional localization error associated with using an 
equivalent magnetic multipole instead of the classical 
current dipole, even if the source is indeed only a cur¬ 
rent dipole. We are motivated by the simpler form of 
the magnetic moment model and the need to reduce 
the search complexity in inverse work. The equiva¬ 
lent CRLB expressions for the current dipole can be 
found in [3]. For the equivalent magnetic dipole, we 
rewrite (8) as 




• m= G(l) • m 


( 11 ) 


The partial derivatives of G w.r.t. location / are given 
by ' 
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where rd is a 3 x 3 tensor product. 


3. Results 


3.1 Performance of ECD vs. EMM for a dipolar 
field 


We first evaluate the performance of the model in 
the simple case of a dipolar field produced by a point 
source. Fig. 1 and Fig. 2 show the CRLB in the posi¬ 
tive x-z-quadrant for the ECD and EMM models 

For both models we assumed ay-oriented dipolar 
source of 10 nA-m intensity, a 35 fT noise standard 
deviation and a 138 axial gradiometer sensor config¬ 
uration located at 10.5 cm from the head origin and 
symmetrically covering the upper hemisphere. The 
contour lines are labeled with the standard deviation 
a of the localization error bounds (in mm). Because 
the EMM has an extra degree of freedom that allows 
the addition of a magnetic dipole, we see a slight 
increase in error in the case of an EMM. The differ¬ 
ence is negligible, however, for shallow regions, indi¬ 
cating that the EMM can safely be used as an 
adequate substitute for the ECD model. 

To check the efficiency of our estimator, we per¬ 
formed Monte-Carlo simulations for both models 
based on a RAP-MUSIC estimator [6]. Several loca¬ 
tions along the z-axis were tested with 100 random 
realizations of noise, and the root-mean-square 
(RMS) error was compared to the CRLB. The results 
in Fig. 3 and Fig. 4 confirm the formulas and demon¬ 
strate the closeness of the CRLBs to the actual RMS 
error results from the simulations. 



Equivalent Current Dipole: CRLB vs. RMS 
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Figure 1: CRLBs for the localization error of an 
ECD in a spherical head model with radially ori¬ 
ented sensors. 
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Figure 2: CRLBs for the localization error of an 
EMM in a spherical head model with radially ori¬ 
ented sensors. 

3.2 Localization error for cortical patch activa¬ 
tion 

To evaluate the performance of the multipolar 
model for distributed sources, three activation 
patches were created on the cortical surface, as shown 
in Fig. 5. Patch #1 comprises 42 dipoles, each 
arranged normal to the local cortical surface, and 
each representing approximately 0.4 mm 2 area of 
cortex. We adjusted the current density in this patch 
to achieve three different source models, equivalent 
respectively to the current dipole, the magnetic 



Figure 3: Monte Carlo RMS error and CRLB along 
the z-axis for the ECD model. 



Figure 4: Monte Carlo RMS error and CRLB along 
the z-axis for the EMM model. 


moment, and the multipole. The magnetic moment 
was formed by explicitly removing the current dipole 
model, then fitting the magnetic moment. Similarly, 
the multipole was formed by explicitly removing the 
magnetic moment, then fitting a multipole. In this 
manner we achieved a “pure” representations of each 
of our models for this patch. For each model, we 
adjusted the source intensity to achieve the same field 
peak of about 100 fT. We then ran 1,000 realizations 
of random noise added to this model, each time local¬ 
izing using the appropriate model in “BrainStorm” 
(http://neuroimage.usc.edu). The RMS errors in 
localizing the true location are shown in Table 1. We 
then repeated this set of models for 100 realizations 
for Patches 2 and 3 in the figure, which comprised 
134 and 137 cortical dipoles respectively, and their 
results are summarized in Table 2 and Table 3 
respectively. 


Table 1: Localization error in terms of CRLB and 
RMS for simulated cortical patch #1. 


Patch #1 

CRLB (mm) 

RMS (mm) 

ECD 

3.85 

3.90 

EMM 

5.24 

5.29 

1 st-order Multipole 

-- 

7.53 



















Figure 5: A tessellated human cortex showing map¬ 
pings of three sources onto the cortical surface. 


Table 2: Localization error in terms of CRLB and 
RMS for simulated cortical patch #2. 


Patch #2 

CRLB (mm) 

RMS (mm) 

ECD 

2.56 

2.53 

EMM 

6.63 

6.67 

1 st-order Multipole 

-- 

14.34 


Table 3: Localization error in terms of CRLB and 
RMS for simulated cortical patch #3. 


Patch #3 

CRLB (mm) 

RMS (mm) 

ECD 

2.73 

2.67 

EMM 

7.40 

7.07 

1 st-order Multipole 

-- 

15.73 


4. Discussion 


Comparing the localization error for the equiva¬ 
lent magnetic moment to that of the equivalent cur¬ 
rent dipole model for a dipolar field, both 
theoretically and through Monte Carlo simulations, 
shows negligible deterioration in localization accu¬ 
racy. It therefore appears safe within the proposed 
multipolar framework to use the EMM instead of the 
ECD to localize true current dipole sources. 

The results for the simulated cortical patches 
show good agreement between the RMS errors and 


the CRLBs for the current dipole and the magnetic 
moment models, indicating the utility of the CRLBs 
in analyzing their performance, without the need for 
extensive Monte Carlos. We note, however, that the 
RMS errors for the multipolar model are in general 
substantially better than the CRLBs would indicate. 
Our multipolar model is unstable in some of its qua- 
drupolar components, which we routinely regularize 
in BrainStorm. The regularization method chosen 
here was truncation of the smallest singular values, 
which introduces a transformation and bias of the 
variables. The CRLBs discussed here are for unbi¬ 
ased analyses only, and our future work will explore 
the more difficult analyses of regularized models. 
Moreover, we expect that further simulations as well 
as human data will reveal the flexibility and effi¬ 
ciency of the multipolar framework in localizing dif¬ 
ferent neural source patterns thereby providing an 
important extension to the traditional ECD model for 
MEG source estimation. 
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